function [rate c f sigmac sigmaf]=run_bounded()

clear all;
load trainingdata;
load testdata;

% Normalizing the feature
data=normalize_feature(data);
testdata=normalize_feature(testdata);

% Shuffle data
order = randperm(size(data,1));
data = data(order,:);

% [beta c f]=sgd_bounded(data,0.1,0);
% sigmac=sigma(c);
% sigmaf=sigma(f);

% % Find best settings
indices = crossvalind('Kfold',size(data,1),5);
options = optimset('MaxIter',100,'PlotFcns',@optimplotfval);
bestsettings=fminsearch(@(settings) crossvalrun_bounded(data,settings,indices),[.1 0.01],...
    options);



[beta c f]=sgd_bounded(data,0.1,0);

sigmac=sigma(c);
sigmaf=sigma(f);

[tp,fn,fp,tn]=test_bounded(testdata,beta,c,f);
rate=(tp+tn)/(tp+tn+fn+fp);
display(['End Processing...']);
